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I. INTRODUCTION 

Amphiphilic molecules in water or oil solutions have 
been intensely investigated over the last decade in a va- 
riety of settings. Amphiphiles in aqueous solution, un- 
der appropriate conditions, have been experimentally ob- 
served to form vesicles with the topology of the sphere, 
torus 1] and even higher genus surfaces [2| ■ Other exper- 
imental studies have focused on the phases of amphiphilic 
films as a function of temperature. R is by now well es- 
tablished that the high temperature fluid phase goes into 
a smectic-C Lpi phase through an intermediate rippled 
P/3' phase. Within the Lp> phase itself there are sev- 
eral other phases characterized by the degree of tilt and 
hexatic bond order 0, This beautiful experimental 
work may provide insight into biological problems such 
as membrane fusion Isf, where it has been argued Q 
that molecular tilt plays an important role. 

A remarkable understanding of the shapes of fluid am- 
phiphilic systems has been provided by physical meth- 
ods based on the Helfrich Hamiltonian 8] and its varia- 
tions 0, . In [ij , the problem of fluctuating smectic- 
C membranes, previously investigated in Il7l| for planar 
films, was addressed and predictions for the shape as a 
function of the elastic constants were presented. R was 
found that toroidal vesicles were favored for some pa- 
rameters. R was assumed, however, that free disclina- 
tions are energetically unfavorable and may therefore be 
ignored unless, as in the case of the sphere, topological 
constraints require them. 

The main result of this paper is that disclinations can 
be energetically favored over a wide range of parameters, 
even when not required by topological constraints. We 
treat vesicles that are topologically torii (closed surfaces 
of genus one). The Gauss-Bonnet theorem for torii re- 
quires a vanishing total disclination charge. We assume 
hexatic order in the tangent plane of the torus, arising 



from an anisotropic liquid phase of molecules with zero 
shear modulus. If {xi,X2) are coordinates on the torus, 
specified in three dimensions by a function R{xi,X2), lo- 
cal hexatic order can be described by a bond angle field 
d{xi, X2) (up to rotations by 27r/6 = 60°) defined relative 
to the local tangent vectors el cx and 62 oc 82 R. On 
the torus, this pair can be chosen to be nonsingular and 
orthonormal everywhere; el ■ ej = 6ij. As discussed, e.g., 
by David [13, the usual hexatic energy on such a curved 
surface can be written 



E = -Ka 
2 



(1) 



where Ka is the hexatic stiffness constant, g^-' {xi,X2) is 
the inverse of the metric tensor 
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and g — det gij . 

The vector-potential-like function U,j[xi,X2) in Eq.JQ) 
describes the geometric frustration which arises when 
vector or tensor order parameters are parallel transported 
on curved surfaces. This "spin connection" can be com- 
puted from covariant derivatives acting on el and €2; the 
"curl" of 17 (when appropriately defined on a curved sur- 
face) is proportional to the local Gaussian curvature |3 ■ 

Disclinations can be inserted into the free energy 
just as in flat space. If TV defects with charges qj = ±1 
are present on the torus, we first define a defect density 
s as a function of a; = (xi, X2), namely 
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As discussed, e.g., in 0,0,1131, minimizing ^ sub- 
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ject to a fixed arrangement of N defects, leads to 



X [six)~K{x)]G{x,y)[s{y)-K{y)] 

K 

+ 2 



Sx^^)H^{x) 



(4) 



where K{x) is tlie Gaussian curvature and G[x, y) is the 
inverse Laplacian on the torus. Up to subtractions which 
eliminate "zero modes" (see Sec. G{x,y) obeys 



< a < Although the {a, 6) coordinate system 

has a clear physical interpretation, most of our results 
are obtained using the conformal coordinates introduced 
in the Appendix [2l| . Upon replacing a by a new angular 
variable f{0<ip< 27r) defined by 



cos (p 



Ri cos a + i?2 



i?2 cos a 



(8) 



one obtains a locally flat metric, which greatly simplifies 
the calculations. 



VlG(x,y) = d,{^g'^d,}G{x,y) 
= S^'Hx-y)/^/^ . 



(5) 



The first term of (0J arises directly from iQJ and repre- 
sents a kind of two-dimensional electrostatics in curved 
space. As discussed below, this electrostatics can lead 
to a lower energy when positive and negative disclina- 
tions are placed at positions on the surface which approx- 
imately match the local Gaussian curvature. To account 
for the bending energy of the surface, we have added the 
second (Hclfrich) term, where the bending rigidity cou- 
pling is K w (1 — 10)kBT for lipid bilayers and H{x) is 
the mean curvature. Defect core energies, which depend 
on short distance physics not accounted for in this con- 
tinuum approach, will be added later. 

Although a representation of the physics of geometri- 
cal frustration like 1^ is possible locally on any smooth 
surface, coordinates for genus one surfaces (the torus) 
can be found which admit such a representation globally. 
One such coordinate system is shown in Fig^ where a 
point on the torus is specified by 



i?(a. 







[Ri + i?2 COS a] cos ( 
-I- i?2 cos a] sin t 
i?2 sin a 



(6) 




FIG. 1: Coordinates (a, 6*) defining the torus. 

Here i?2 < -Ri so that the torus does not intersect itself. 
The Gaussian curvature associated with © is a function 
of a only. 
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(7) 



Note that K is positive on the outside wall of the 
torus (— ^ < a < -1) and negative on the inside wall 
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FIG. 2: Patch of positive Gaussian curvature on the torus 
(shaded) subtending an azimuthal angle A6. The surface of 
the torus is divided into regions of positive and negative Gaus- 
sian curvature by the circles labelled A and B, on which the 
Gaussian curvature vanishes (after D. Hilbert and S. Cohn- 
Vossen, Geometry and the Imagination (Chelsea, New York, 
1952). 

The "electrostatic" energy term in Eq.Q favors ap- 
proximately charge neutral configurations, with discrete 
disclination charges cancelling the smeared out Gaussian 
"curvature charge". Although the full calculation (with 
core energies taken into account) is subtle (see Sec lIVII . 
it is interesting to estimate how many additional discli- 
nations might be accommodated on a torus. In the torus 
shown in Fig|21 the solid and dashed circles divide the 
surface into regions of positive and negative Gaussian 
curvature. Consider a wedge of angular width on 
the outside wall of positive Gaussian curvature. The net 
"curvature charge" associated with this region is 



AK = del da^K{a) ^2{M) , (9) 

Jo J-tt/2 



where we have used Eq.Q and ^ — 
R1R2 + ^ cos . Upon equating AK to 27r/6, 
the charge of a single disclination, we find that 
AO = 27r/12, independently of i?i and i?2- Thus 
2it/ A9 = 12 positive disclinations would be required to 
completely compensate the curvature of the outer wall. 
Similarly, 12 negative disclinations would be required 
to completely compensate the negative curvature of the 
inner wall. This simple argument neglects core energies 
and interactions between disclinations, effects which will 
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cause the preferred number of defect pairs to be less 
than 12. 

The thermodynamic limit for torus dimensions Ri and 
i?2 corresponds to the limits Ri/ao oo and i?2/ao ^ 
oo, with the aspect ratio r = R1/R2 fixed, where oq 
is a microscopic length scale such as the particle spac- 
ing. Upon optimizing r with a defect-free hexatic tex- 
ture on the torus, one finds that a nonzero hexatic stiff- 
ness constant Ka pulls r above the Clifford torus value 
r = \/2 appropriate for liquid torii 1 so that the result- 
ing shape looks more like a bicycle tire (see Sec. lIVt A'). 
Our results for defect energies are presented for fixed r 
and a given number oiM ^ R^R^Ia^^ of microscopic de- 
grees of freedom. It would be straightforward, however, 
to use the methods described here to optimize over both 
r and possible defect configurations for fixed M . 

Although we find that disclinations are always unfavor- 
able in the "thermodynamic limit" of large M, the critical 
value M = Mc required to suppress them in the ground 
state is surprisingly large. Indeed, for r — \/2, this num- 
ber is of order VS^^\ (see Fig ll2|) . As r approaches 1 from 
above, i.e., in the limit of "fat" torii, we find that Afc 
exhibits a remarkably strong divergence [see Ea. H49|l ]. 



as treated in Sec. llVt C. Because is so large, the 
toroidal vesicles of Ref. would be quite likely to have 
disclinations present in the ground state if hexatic or- 
der were present. Indeed, for R\ — 5/i (roughly the size 
of a red blood cell), Ri = rRi « 1.4i?i and ao = 20 A 
(typical lipid spacing in vesicles), we have (see Sec lIVt Bl 
M = « 4 X 10* which is much less than the crit- 

ical value Mc ~ 10^^ required to suppress disclinations. 
As discussed above, the interaction between hexatic or- 
der and the Gaussian curvature of toroidal vesicles leads 
to r > V2 and a smaller value of Mc- Vesicles with lipids 
in a liquid state could provide, however, a toroidal tem- 
plate with r = \f2 for adsorbed colloidal particles, similar 
to the spherical "coUoidosomes" studied by Dinsmore et. 
al. • It may be possible to use polymerizable lipids to 
permanently fix the template aspect ratio at r = \/2- It 
would be quite interesting to study (with, say, confocal 
microscopy) both hexatic and crystalline order in col- 
loidal particle arrays jlS^ adsorbed on such a template, 
as has already been done for colloids on spherical water 
droplets in oil 0| . The colloid radius would play the role 
of a microscopic scale ao ^ i?i , R2 in this case. Discli- 
nation defects in a crystalline ground state might well be 
accompanied by grain boundaries |l9j| . 

Although we focus here on hexatic order in toroidal ge- 
ometries, similar results should apply to XY-like models, 
as would be appropriate for vesicles composed of lipid 
bilayers with tilted molecules 0, . Our results are rel- 
evant as well to two-fold nematic order on a toroidal 
template. In both cases, we expect qualitatively similar 
phenomenon, such as defects in the ground state, unless 



the total number of degrees of freedom exceeds a critical 
value. A precise equivalence is possible in the one Frank 
constant approximation As discussed in Sec. IIVI 

defects in the ground state are more likely for fat torii 
in the case of nematic (and hexatic) order. Interesting 
results related to those here have recently been obtained 
for "corrugated" topographies, which are flat at infin- 
ity and for which there is also no topological necessity 
for defects in the ground state |3- Specifically, it has 
been shown that defect pairs lower the energy of a hex- 
atic phase draped over a Gaussian "bump" for sufficiently 
large ratio of height to width. In this case, defects remain 
an important feature of the ground state even when the 
number of degrees of freedom M tends to infinity. 

The organization of the paper is as follows: an analogy 
to an electrostatic problem, aimed at providing a more in- 
tuitive picture of the physics of curvature-induced defect 
unbinding, is introduced in Sec^l In Sec. Illll the inter- 
action among defects is worked out in detail for toroidal 
topology. Predictions for the total number of defects are 
provided in Sec. IIVI The effects of both temperature and 
shape fiuctuations are discussed in Sec. 

II. ELECTROSTATIC ANALOGY 




FIG. 3: A plus (filled circle) and a minus (filled square) pair 
of charges between the plates of a two-dimensional circular 
capacitor is analogous to a plus-minus pair of disclinations on 
a torus with the identifications given. 

It is useful to illustrate the physics of defects on the 
torus with a simple electrostatic analogy, illustrated in 
Fig- El A positive and negative charge are placed between 
the plates of a capacitor with circular cross section. The 
plus-minus pair are the analogs of a plus-minus disclina- 
tion dipole initially located on a zero curvature line of 
the torus (see Fig. 0J). We assume here that the charges 
are extended over a core radius ao, which plays the role 
of a minimum separation. We expect that this core ra- 
dius is related to the mean particle separation on the 
torus. The net charge -l-Q and —Q on the capacitor 
plates represents the Gaussian curvature in Eq.Q), in- 
tegrated over the regions of the torus where it is positive 
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and negative, respectively. The linear charge density on 
the plates of the capacitor is thus given hy p± — Q / L±, 
where L± — 27r(i?i=Fi?2). The lengths Ri and R2 for this 
capacitor correspond to the torus radii defined in Fig. ^ 



Eq.|T2l 

ance of 
ing d = 



illustrates very clearly the appear- 
two preferred locations; these be- 
ao (where £ w 2£'c) and d ^ R2 



Ri +B.2 

R1-R2 



The 



[where £ = 2Ec + ln(i?2/ao) - eQ In 
relative strength of the first two terms will determine 
the preferred location of the charge dipole. By taking oo 
small enough, for Ri > R2, the attractive charge-charge 
term will dominate and defects will not be favored. For 
any finite ao, on the other hand, defect unbinding will 
be favored in the limit R2 — > R^ ■ This is precisely what 
we find, up to numerical constants, in our treatment of 
disclination unbinding on the torus. 



FIG. 4: Illustration of the calculation discussed in the text: 
a plus (filled circle) and a minus (filled square) form a discli- 
nation dipole on one of the two circles of zero Gaussian cur- 
vature. They are then pulled apart until they reach the max- 
imum curvature line (plus) and minimum curvature line (mi- 
nus). 



The competition between the work done by the electric 
field of the capacitor (representing regions of +/— Gaus- 
sian curvature on the torus) and the attraction of op- 
posite sign charges will dictate whether the disclinations 
separate or remain tightly bound at separation a. Since 
the energy of a configuration with excess bound charges 
exceeds that for no excess charges by an amount of or- 
der two core energies, this criterion determines whether a 
separated plus-minus pair is present in the ground state. 

In this two-dimensional geometry, the electrostatic en- 
ergy of the two charges with separation 2c? and charge e 
is given by 



£1 ^ ln(d/ao) + 2Ec 



(11) 



where Ec represents a self-energy of an isolated charge, 
corresponding to the core energy of a disclination. There 
is an additional electrostatic force pulling the charges to 
the capacitor plates, which leads to an additional energy 



£2 = Qe In 



Ri-d 



Ri +d^ 
The total energy is then 

£ = 2TTe^Q-QeC + 2Ec 

where 

Q = ^ Hd/ao) , 



(12) 



(13) 



C = In 



Ri+d 
Ri-d 



(14) 



The functions Q and C have a similar form to those we 
find in the exact calculation on the torus (see Ea. (|31|l ). 



III. INTERACTING DEFECTS ON CURVED 
SURFACES 

A. Green's function on the Torus 

For arbitrary coordinates x = (xi, X2) on the torus, the 
key object in the energy Eq. Q is the inverse Laplacian, 
which is the solution G(x,x') to the equation 



AG(x,x') =<5(x,x') 



1 

A 



(15) 



where the Laplacian is defined as in Eq.(|SJ, (5(x, x') 
--^(5(x — x'), and A is the area of the torus. 



A = y d'^Xy/^ = in'^RiR2 



(16) 



A constant is subtracted from the (5-function to eliminate 
a "zero mode" which changes the area of the torus. As 
shown in the Appendix, the metric of a torus, with a 
suitable change of coordinates, is flat modulo an overall 
conformal factor [see Ea. (|65|l ]. The inverse Laplacian 
may be computed by considering the flat metric (where 
the conformal factor is identically one) and its associated 
inverse Laplacian G(x, x'), from the formula 



G(x,x') =G(x,x') 



A 

A I A 



,My) G(x,y)+G(y,x') 



V5(y)V5(?)G(y,y')(17) 



As can be checked straightforwardly, G(x, x') solves 
Eq. I|15|) as well as satisfying 



J d2xy^G(x,x') = 
J d2xV5(x')G(x,x') = 



(18) 



Thus Ea. H17|) is indeed the inverse Laplacian, where the 
conditions (|18|l ensure overall "charge neutrality" for any 
disclinations present on the torus. 
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The coordinates for the torus are shown in Fig. ^ 
Upon making the change of variables a — *■ (/3 via 



cos a 



Ri cos — 



D — D , described in the Appendix, the in- 

/tl — it2 cos ' ' 

verse Laplacian G(x, x') in conformal coordinates can be 
computed by first solving 



— (sinhp9g 



1 



sinh p 



dl)G{e,v\e',^') 



(27r) 



(19) 



where sinh p = \/r^ — 1 and r is the aspect ratio i?i /R2 
of the torus. The straight-forward solution 



47r 



■In 



(sinh p){9 



G(")(0,^|e',^') 



sinh p 



(20) 



satisfies AG'^°) ~ 5{9 — 9' ,tp — (p'), but is not periodic, 
i.e., 

G^"\9 + 27Tk,^\9',^') ^ G'^°'>{9,(p + 2TTn\9',ip') 

^ G^'H9,^\9',^') , (21) 

for arbitrary integers k and n. In addition the Laplacian 
acting on G'-^\9, (p\9' , ip') fails to give the constant term 
on the right hand side of Eq. H19(l . Both these deficiencies 
are remedied by defining 

G{9,ip\9' ,p') = --J-V ln[sinhp(6l + 27rfc)2 
27r ^ — ' 

k,n 

(^ + 27rn)2] +C(p) , (22) 



sinh p 



where the constant C{p) is determined by imposing that 
the inverse Laplacian exhibits the correct short-distance 
singularity (following from Ea. (|20|l 'l in the limit 9^9' 
and kp ^ if)' . Standard analytical techniques [2^ l23 | 
allow one to perform the sum indicated in Eq. (|22|l . giving 



G(0,^|0',^') = -—In 



1 



sinh p\d\(^ 



2tt 



2 sinh p \ 2n 



4:n^\'i]{T) 
2 



(23) 



where r = . ' . The functions i?i and 77 are the Theta 

Sinn p ^ ' 

function and Dedekind i] function, respectively, defined 
by HI 



n— +00 



and 



/2'f 2TTw{n-l/2) 



n— — OQ 



„2irir/24 



(24) 



(25) 



The inverse Laplacian of Eq. H17|l thus becomes 
G{9,p\9',p')^G{9,p\9',p') - 

y ^ — —[cos{np) + cos{nip')] 



{2'kY \ sinhp 

\ n—l 



+ 



cosh 



- — [cos(ni^) 4- cos(ni^')] + 

shp ^-^ n I 

n—l / 

2 f 1 ^ 



-2np 



(27r)2 \^sinhp^^ 



2 - e^2"^tanhp - 



, (26) 



cosh p ^-^ n cosh p ^-^ 

n—l m—1 

which can then be used to evaluate Eq. Q). 

B. Energetics of Defects on a Torus 

The total energy Eq. also contains a bending rigidity 
term which, for a torus with aspect ratio r — R1/R2, is 



(27) 



If we were to minimize this term alone, we would find 
r = v^, the so-called Clifford torus. The Gauss-Bonnet 
theorem for a torus reads 



(28) 



and, with our choice of Green's functions, Eq.(^) insures 
that the sum of disclination "charges" qi satisfies 



N 



(29) 



so that, as previously noted, no defects are required topo- 
logically. The defect charges here take on the values 
qi = ±1,±2, . . ., with qi = ±1 corresponding to the 
elementary defects with rotations ±27r/6 in the hexatic 
order parameter. 

With the Green's function in hand, the hexatic energy 
of a set of disclination charges on a torus (the first part 
of Eq.(E3)) can thus be written explicitly: 



(30) 



where the defects interact with each other according to 

^^^^ ^_ iiVi~Vj) i \\2^ 

S(x,,Xj) 



-L , i \"1\ 2tt ^ 27r sinhp ' sinhp''! 



+ 



1 



2 sinhp \ 2n 



Pi - pj 



I ' ^ smh p ' I 

2 



(31) 
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and interact with the background Gaussian curvature 
"charge" according to 



/:(x) = In 



1 



^ cosh p — cos 

The "spin wave" part of the frustrated hexatic energy 



(32) 



(33) 



is present even without defects and supplements the 
bending rigidity term Eq. I|27|) . The core energy term 
in Eq. (|30|l will be considered in more detail in Sec. IIVI 



C. Energetics of Defects on a Sphere 

It is instructive to compare our results for toroidal vesi- 
cles with the corresponding results for spherical vesicles. 
The interaction potential for a sphere is 



1 — cos j3i. 



(34) 



where /3y is the geodesic distance, for a sphere of unit 
radius, between points x^ and Xj. The function £(x) and 
the constant V can both be set to zero on the sphere. 

The Gauss-Bonnet theorem for spherical topology 
reads 



yielding the constraint 



N 



12 



(35) 



(36) 



IV. GROUND STATES OF HEXATIC 
TOROIDAL VESICLES 

Before presenting a detailed analysis of the implica- 
tions of the Hamiltonian (|30l) for defects on a torus we 
must confront the issue of core energies. The short- 
distance behavior of the defect potential Q (Ea. l^'Sl\i ) im- 
plies that a plus-minus pair of disclinations located as in 
Fig0|on a circle of fixed azimuthal angle 9 and geodesic 
separation 2d have an energy 



E : 
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liiiR2/ao) + Vi{d/R2)]+2E, 



(37) 



where, after absorbing various constants, V7(d/i?2) may 
be viewed as the interaction potential between the two 
disclinations and the term Ec (reflecting short-distance 
physics) has been added by hand. For <C 2d <C i?2, 
Vi{d/R2) ~ ln(2d/i?2)- The torus radius i?2 therefore 
drops out in this limit and we recover the result E = 



-fCATr/lS ln(2d/ao) -I- 2Ec for a disclination pair in flat 
space. The disclination core radius ao is the distance (of 
order the spacing between molecules) at which continuum 
elasticity breaks down. We will henceforth assume that 



Ec ~ cKa , 



(38) 



where c is a dimensionless parameter characterizing the 
size of the disclination core energy. In our numerical 
calculations we take c = 0.1. 

If the torus is coated with M particles, an effective 
average particle spacing ap (assuming a local triangular 
lattice for simplicity) may be defined from the area per 
particle 



V3 



M 



Att^RI r 
M 



(39) 



We shall assume that oq ~ ap, so Ea. (|39|l relates the 
total number of particles M to the minimum plus-minus 
pair separation appearing in Ea. H37|) . 

Upon including both the bending rigidity and the core 
energies, the total energy for N charge ±1 defects on a 
rigid, undeformed torus takes the form 



IT 10^ 



N 



i=l j>i 

47r2 



1 



2TT^r^ , , 
K , „ (40) 



E depends on the lattice spacing oq through the core 
energy Eq. as well as the constraint that the defect 
spacing cannot be smaller than ap. In our calculations 
the latter constraint is accounted for by forbidding any 
two defects from approaching within the distance ap or, 
alternatively, by assuming that any two defects closer 
than ap merge into a single defect with total charge the 
sum of the individual charges of the two defects. 

Only the potential (Q) and curvature (£) terms de- 
pend explicitly on the defect positions and charges. For 
fixed R\ and i?2, therefore, the total number of defects 
in the ground state is independent of both the bending 
rigidity and the constant T) in Ea. H33|l . 

The defect-defect interaction is determined through 
the Q-function defined in Eq. (|31|l . For two opposite sign 
defects the energy is attractive for all separations at con- 
stant 0, as shown in Fig.[S| If only this term were present, 
the attraction would bring both charges as close as pos- 
sible, binding all disclinations into dipoles which have a 
higher energy than a defect-free configuration. Thus if 
no other terms were present, the ground state would be 
defect free. 

The defect-curvature interaction L favors the appear- 
ance of additional defects. This term acts like an elec- 
tric field pulling the positive/negative disclinations into 
regions of positive/negative Gaussian curvature, respec- 
tively, similar to the electrostatic analogy discussed in 
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FIG. 5: Various contributions to the energy (in units of Ka/2) 
with aspect ratio r = y/2 for a disclination dipole separated 
an angular distance a along the path shown in Fig. |1] The 
long dashed line is the energy in the absence of defects. The 
dashed-dotted line is the defect-defect interaction as a func- 
tion of the separation of the charges. The dashed line is the 
curvature-defect interaction energy as a function of separa- 
tion, and the continuous line is the total energy. The core 
energy contribution, computed from Eq. 1381 is very small on 
this scale. The bending energy is subtracted for clarity. 

Sec^ As shown in Fig. [SJ if this were the only term 
present, isolated plus and minus disclination charges are 
always energetically favored, with the lowest energy aris- 
ing when they are located at the regions of absolute max- 
imum Gaussian curvature. 

The total energy is a competition between defect at- 
traction and curvature-induced unbinding, as shown in 
Fig. El where the potential shows a double well, corre- 
sponding to the attractive defect-defect interaction min- 
imum and the attractive Gaussian curvature minimum. 
In Fig. the curvature field dominates and additional 
defects appear in the ground state, but the general result 
for whether defects should or should not be expected is a 
function of the core energy Ec, the hexatic stiffness Ka, 
the torus aspect ratio r — R1/R2, and the ratio i?2/ap 
of macroscopic to microscopic cutoff. 



A. Defect-free hexatic toroidal vesicles 



which read in the limit of large and small hexatic stiffness 

El, 

r = V2 , Ka « K 

r = , Ka» >i . (43) 

V 2k 

If the hexatic stiffness is much smaller than the bend- 
ing rigidity, the Clifford torus (r — ^/2) is the optimal 
geometry. If, on the other hand, the hexatic stiffness 
dominates, then a thin torus, similar to a bicycle tire, is 
optimal. This picture changes when defects are included. 



B. Ground states of defective hexatic toroidal 
vesicles 

The general ground state for arbitrary aspect ratio may 
be determined numerically using Eq. H4U|I . For simplic- 
ity, we compare the energies of configurations with and 
without defects for a fixed aspect ratio r. 

We first performed the following calculation: a set of 
N unit charge disclinations {N/2 positive and N/2 nega- 
tive) are placed in opposite-sign pairs on a circle of zero 
Gaussian curvature. Each pair is then pulled apart at 
constant azimuthal angle 9 until the plus (minus) discli- 
nations reach the outer (inner) rim of the torus respec- 
tively. This is illustrated in Fig. 0] for the simplest case 
N ^2. 




FIG. 6: The total energy (in units of Ka/2) for aspect ratio 
r = -^5/4 for varying numbers of defects and ap — -R2/4. 
The bending energy at fixed r is subtracted off. The discli- 
nation core energy is taken to be O.IK a, which is 0.2 in the 
above units (corresponding to c = 0.1). 



In the absence of defects the total energy following 
from Eq. (j^Ol is 



2tt^Ka 27r2r2 

E = -j^=^= + K = 



(41) 



a result first obtained in The optimal value of r 

resulting from minimizing this energy is the solution of 

(r2 - 1)3/2 _ ^(^2 _ -^^ ^ —r(r^ - 2) = 0, (42) 

Ka 



We discuss our results for two regimes of disclination 
core energy. Core energies corresponding to c-coefRcients 
less than 1/10 have almost no effect on the energy bal- 
ance since the elastic energy is already quite large. For 
definiteness the small core energy regime will be illus- 
trated for c = jij. Our results are shown in Figs. |H1 E| 

and|Slfor the three aspect ratios: r = ■\/5/4, r ~ \/2 (the 
Clifford torus) and r = 2.69. We have set i?2/ap = 4, 
so that the number of degrees of freedom (see Eq. (|39|) ) 

on these torii is M = ^rl^ f ~ 800, 1000 and 2000 
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FIG. 7: The total energy, as in Fig. El but for the ChflFord 
torus which has aspect ratio r = y/2. 



respectively. In each case the addition of defects lowers 
the energy. Note that the energy at maximum separation 
{a = tt) first decreases, and then increases with N. The 
optimal number of defect pairs {N/2) is 5, 6 and 7 for 
Figs. El and |H1 respectively. This number is less than 
the naive estimate of 12 in the Introduction due to re- 
pulsive defect interaction energies on the inner and outer 
walls of the torus. 

The typical situation for large core energies is illus- 
trated for c = 1. It is only for core energies of this 
order that we find qualitatively different behavior from 
the small core energy regime. Our results are plotted 
in Fig. Mir = ^^S- ^ = be- 

cause i?2/ap = — — 32 in these plots, we now have 
M « 52, 000 and M « 67, 000 respectively. Even with 
such a large core energy, defects are present in the ground 
state with a preferred number of pairs N* /2 ~ 3 for 
Fig. Inland a smaller number for Fig.^1 With this value 
of i?2/ap the torus is defect free for larger r. 

Note that for fixed aspect ratio the preferred number 
of defects drops for larger numbers of particles when the 
core energy is large: compare Fig. |51to Fig. [71 To study 
this point further we have determined the total num- 
ber of defects in the ground state as a function of the 
number of particles. Our results are shown in Fig. ^2 
using M — ^r(i?2/ap)^ and assuming oq = ap. Al- 
though a torus always becomes defect free in the ther- 
modynamic limit i?2/flp ^ oo (with r = R1/R2 fixed), 







r=2.692S8 

N=a 






N=24 






^~ — — N=4 






^ N=20 






" - N=8 






i:::;;-^^ n=i6 










FIG. 9: The total energy (in units of /('a/2) for a torus with 
aspect ratio r = 1^5/4 for varying numbers of defects with 
R2/ap = IOO/tt. The disclination core energy is taken to be 
Ka which is 2 in the above units (corresponding to c = 1). 
The bending energy at fixed r is subtracted off. 



torii with moderate aspect ratio only become defect-free 
for a remarkably large number of particles, which may be 
as large as 10"'^^ for r = ^/2\ 

To make the last point more transparent, the criti- 
cal number of particles M^, above which defects are no 
longer favored, is plotted as a function of the aspect ra- 
tio in Fig. ^1 As the aspect ratio r = i?i/i?2 ^ l"*", Mc 
diverges, suggesting that any torus will possess defects if 
sufficiently fat. We provide an rough analytic argument 
along this line in the next subsection. 

Although fat torii (with r > 1) tend to favor defects 
in the ground state, the maximum number of defects fa- 
vored for a given aspect ratio is a subtle question. To see 
this, note that the constrained minimization discussed 
above leads to a ring of iV/2 positive disclination charges 
on the outer wall of the torus, and a smaller ring of Af/2 
negative charges on the inner wall. As r — s- 1+, negative 
charges end up being very close in the final configuration, 
with a considerable repulsive energy cost. 

To allow the system to reduce this energy, we have 
considered the modified calculation illustrated in Fig. IT!^ 
The initial configuration starts from both circles of zero 
curvature, alternating plus-minus defect pairs between 
them. When the separation variable a = tt the config- 
uration is the same as in the previous case. A typical 




FIG. 8: The total energy, as in Fig. |S1 but for aspect ratio 
r = 2.6926. 



FIG. 10: The total energy, as in Fig. |2l but for aspect ratio 
r = V2. 
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M (Number oi particles) 



FIG. 11: The preferred number of defects as a function of 
the total number of particles for three aspect ratios. The 
disclination core energy is taken to be O.IKa- The dashed 
line corresponds to the configuration in Fig. 0] 




FIG. 13: Illustration of the calculation discussed in the text: a 
plus disclination (filled circle) and a minus disclination (filled 
square) form a defect dipole on one of the two zero-curvature 
circles of the torus. They are then pulled apart in the direc- 
tion of the maximum curvature line (plus) and the minimum 
curvature line (minus), respectively. 



plot of the energy is shown in Fig. In this plot, the 
minimal configurations (with a ^ n) correspond to buck- 
led or staggered rings of defects, displaced from the cir- 
cles of maximum or minimum curvature, as illustrated in 
Fig. El The final position is a compromise between the 
"electrostatic" repulsion and the attraction to the Gaus- 
sian curvature basins. With the extra degrees of freedom 
allowed by buckling, the energies are always lower than 
the simple ring configurations of Fig. Eland Fig. [51 Note 
that staggering allows more defects to be squeezed into 
the ground state: the optimal number N* /2 of defect 
pairs is 7 in Fig.^^ as opposed to N/2 « 4 — 5 in Fig. El 





r (aspect ratio) 



FIG. 12: The critical number of particles, above which defects 
are no longer favored, as a function of the toroidal aspect ratio 
for vanishing core energy (blue) and c = 0.1 (green). The 
analytic estimate for c = 0.1 is plotted as solid red squares. 



FIG. 14: Plot of the energy for the path described by Fig. 1131 
The disclination core energy is taken to be O.IKa- The aspect 
ratio is r = -^5/4 and R^/a — 50. The optimal number of de- 
fect pairs in the final "buckled ring" configuration is N* /2 — 7 
for this M « 127, 000 particle configuration. 



C. Analytical argument for defects in the ground 
state for aspect ratio r near 1 

Let us consider a +/— disclination dipole on a zero 
curvature circle of the torus. Imagine slowly pulling the 
dipole apart until the individual +(— ) disclinations reach 
the outer(-|-)/inner(— ) rim of the torus respectively. The 
total energy in this configuration is dominated by the 
attraction of each defect to the corresponding region of 
same sign curvature since the defects are too far apart for 
the defect-defect interaction to be important. The total 
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FIG. 15: The typical ground state configuration for parame- 
ters that favor defect proliferation. The arrows indicate the 
displacement of the equilibrium defect position from the max- 
imal curvature circles. 



energy following from Ea. (|30|l is therefore 



(44) 



where we have set ipi — and (^2 = in Ea. (|32|) . Upon 
approximating the defect-pair energy by its flat space 
value 



£fir] ~ — Ka In I — - 
18 \ao 



(45) 



- the agreement is excellent. This result also establishes 
that excess defects are present in the ground state for any 
fixed particle number provided the torus is sufficiently 
fat. Defects could thus be an important feature of hexatic 
textures for realistic vesicle sizes. 

It is interesting to generalize these formulae to p-fold 
symmetric order parameters |ll| on the surface of a torus. 
Here, hexatic order corresponds to p = 6, nematic order 
to p = 2 and tilt order top = 1. A hypothetical "tetradic 
phase" with a four-fold liquid crystalline symmetry psl ] 
would correspond top — A. The generalization of Ea. H46|l 
for a minimally charged defect-antidefect pair with p-fold 
symmetry reads 



2tt 

£ = KaIyi 

P 



1 



1 



2tt 



R 



+ ^KAln[ — ]+2E,. (50) 



ao 



This yields a critical particle number 



Mr 



87r2 
V3 



exp 



ttKa 



r - 1 



2p 



(51) 



The critical number of particles above which defects no 
longer appear in the ground state is therefore lower for 
coatings of the torus by textures of lower symmetry. 
Since Ri must exceed R2 by an amount of order ao for a 
physical torus (see Fig.P), (r — l),„in = (00/^2)- Hence, 
Afc diverges like (i?2/ao)^^ in the limit of an extremely 
fat torus. Upon noting that M ^ (^2/10)^1 we see that 
typically M <C Mc whenever i?2/ao ^ 1. Thus, defects 
are an inevitable part of the ground state for sufficiently 
fat torii in all cases, except possibly for p = 1. 



we find a total energy 



E--\KAln{'-^ 
6 V r — 1 



18 V ao 



2Ec . (46) 



If we assume, consistent with our numerical evaluation of 
Q, that the constant correction to Eq.|03J) is negligible, 
then Ec can be interpreted as a microscopic disclination 
core energy. Eq. (|4(j|l changes sign for 



ao 



exp 



-mEr 

ttKa 

\2 



r + 1 
r - 1 



(47) 



Using M — (i?2/ao) , we conclude that defects are 
favored for 



V. TEMPERATURE AND SHAPE 
FLUCTUATIONS 

A. Connection with Two Dimensional Melting 

The renormalized Frank constant for a film in the hex- 
atic phase has the temperature dependence |25| . 



Ka{T) e+jT) 

knT a^ 



(52) 



where ^+ is the correlation length. The correlation length 
itself behaves in the neighborhood of the hexatic to fluid 
transition temperature T/ like 



M < Mr 



Stt^ r -72Er 



1 



(48) 



For the representative value Ec — Q.IKa, we therefore 
find 



Mr « 4.6r 



1 



r- 1 



12 



(49) 



A comparison with our numerical results for Mc for both 
vanishing core energy and Ec = O.IK a is shown in Fig ll2l 



C4 



exp ■ 



(53) 



The bending rigidity has been shown to have a much 
weaker temperature dependence |23, ll^]. Near the 
hexatic-fluid transition, therefore, the ratio of Ka/k di- 
verges, which should produce larger values of r. 

For toroidal vesicles these results change in two impor- 
tant ways: both the finite size and the Gaussian curva- 
ture of the torus must be taken into account. The fi- 
nite area of the torus limits the growth of the correlation 
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length, viz. 



VII. APPENDIX 



or equivalently, 



Ka{T) 



< 



(54) 



(55) 



so that the Frank constant no longer diverges. It is possi- 
ble that the effects of Gaussian curvature will even limit 
Ka/k to smaller values. As discussed in the Introduc- 
tion, it may be possible to "freeze-in" an aspect ratio 
r « -^2 by using lipid bilayers with only short range or- 
der as a toroidal template. 



B. Fluctuating Hexatic Membranes 

In [2^ the properties of a fluctuating hexatic mem- 
brane were investigated. It was found that the long- 
distance behavior is governed by a new fixed point, char- 
acteristic of a crinkled phase intermediate between a 
crumpled and a rigid phase. Within a large d expansion, 
the new fixed point has the property 



Ka 

K 



Ad 

y 



= 4 at d = 3 



(56) 



It can be shown that for the value of r corresponding to 
this ratio of elastic constants, additional defects should 
be present. The aspect ratio as a function of the elas- 
tic constants for a defect free configuration gives (see 
Eq.gnil) 



V2, 



(57) 



a Clifford torus, which we have shown contains additional 
defects in the ground state. 
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In the angular coordinates {6, a} of Fig. (with < 
6 < 2tt, < a < 2tt) the parametrization 



X — {Ri + i?2 cos a) cos 9 , 
y = (i?i -I- i?2 cos a) sin 9 , 

z — R2 sin a 



(58) 



defines a torus as the locus of points (x, y, z) that satisfy 
{^■s/x"^ + 2/2 — Rij +z'^ — R2^. The dimensionless aspect 



ratio r 



R2 



(59) 



is constrained to be greater than one for torii which do 
not self- intersect. The metric is given by 



ds2 ^ + cos afdO'^ + da^ } 



(60) 



Upon introducing a new angle variable f (0 < ip < 2tt) 
via 



cos a 



r cos If — 1 
r — cos (p 



(61) 



Eq.(jnH|) becomes 



y 



a sinh p cos 9 
cosh p — cos ip 

a sinh p sin 9 
cosh p — cos (fi 
a simp 



cosh p — cos ip ' 
where a and p are defined by i?2 = a/ sinh p and 

r = coshp (^p = In |r + ^/r^ — l|^ . 

In these coordinates, the metric Ea. (|50|l becomes 

— 1 \ ^ / .„o d(p2 



(62) 
(63) 



(64) 



ds^ = Ri 



r — cos p 



d9' 



1 



(65) 



The metric is now conformally flat, i. e. up to a Lp- 
dependent multiplicative prefactor this is the metric of 
a plane with (rectangular) periodic boundary conditions. 
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